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Abstract 

We show how the hnearized equations of motion of any dissipative current are determined by the 
analytical structure of the associated retarded Green's function. If the singularity of the Green's 
function, which is nearest to the origin in the complex-frequency plane, is a simple pole on the 
imaginary frequency axis, the linearized equations of motion can be reduced to relaxation-type 
equations for the dissipative currents. The value of the relaxation time is given by the inverse of 
this pole. We prove that, if the relaxation time is sent to zero, or equivalently, the pole to infinity, 
the dissipative currents approach the values given by the standard gradient expansion. 
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I. INTRODUCTION 



During the last decade, relativistic fluid dynamics has been used as one of the main 
theoretical tools jl| in the study of the hot and dense matter created in heavy-ion collisions 
at RHIC and, most recently, at LHC While the properties of non-relativistic fluids are 



now well understood 



3[, the non-trivial physical consequences imposed by causality and 
stability of relativistic fluids are still under intense theoretical investigation 0, Q . 

In non-relativistic Navier-Stokes theory, the dissipative currents, such as the bulk viscous 
pressure, the shear stress tensor, and the heat flow, are assumed to be linearly proportional 
to the fluid-dynamical forces, such as spatial gradients of fluid velocity, temperature, and 
chemical potential. The constants of proportionality are the bulk viscosity, the shear vis- 
cosity, and the heat conductivity. Navier-Stokes theory can be extended by considering 
higher-order gradients of fluid velocity, temperature, or chemical potential, leading to the 
Burnett equations (including second-order gradients), the super-Burnett equations (includ- 
ing third-order gradients) etc. jsl. A systematic derivation of these equations is provided by 
the gradient expansion. The first-order truncation of the gradient expansion, i.e., Navier- 
Stokes theory, is stable, but not causal, as it allows for propagation of signals with infinite 
speed 4|. Higher-order truncations suffer from the Bobylev instability [7], and are thus 
neither stable nor causal. 

Eckart, and later. Landau and Lifshitz were the first to attempt a relativistic formulation 
of fluid dynamics p|, O]. Their derivation was based on a relativistically covariant exten- 
sion of traditional Navier-Stokes theory. The resulting equations coincide with a first-order 
truncation of a relativistic formulation of the gradient expansion. Unlike the non-relativistic 
case, however, already the first-order truncation of the gradient expansion, i.e., the rela- 
tivistic generalization of Navier-Stokes theory, is not only acausal but, unfortunately, also 
unstable 



Israel and Stewart jo] were among the first to formulate a relativistic theory of fluid 
dynamics that respected causality and was potentially stable. Their formulation was a rela- 



tivistic extension of Grad's theory [10|, where the fluid-dynamical dissipative currents appear 
as dynamical variables which relax to the values of Navier-Stokes theory on characteristic 
time scales, usually referred to as relaxation times. Thus, unlike for the gradient expansion 
the dissipative currents in Israel-Stewart theory (IS) do not have to be zero in the absence 



of gradients. Instead, they decay to zero on the time scales given by the relaxation times. 

One of the features of Israel-Stewart theory is that it reduces to Navier-Stokes theory in 
the limit of vanishing relaxation times. In other words, in Navier-Stokes theory the dissipa- 
tive currents relax instantaneously to the values given by the fluid- dynamical forces, which 
leads to a violation of causality. For the theory to be causal it is therefore necessary that 

;he relaxation times assume a non-zero value, but this is not sufficient. It was shown in Ref. 

sl for the case of bulk and shear viscosity that causality imposes a stronger constraint for 
IS theories: the ratio of the relaxation times to the viscosity coefficients must exceed certain 
values j^. For instance for the case of shear viscosity only, the ratio of shear relaxation time 
Tt, to shear viscosity 77 must obey the relation 

4 

> TTT^-^ , (1) 



+ P) - 3(1 -C2) 

where e, P, and Cg are the energy density, the thermodynamic pressure, and the velocity of 
sound, respectively. It was also shown in Ref. js] that, for relativistic fluids, the causality con- 
dition (II]) also implies stability of the fluid-dynamical equations. Any physically meaningful 
theory of fluid dynamics must not be intrinsically unstable, and any relativistic generaliza- 
tion of such a physically meaningful theory must be causal. Therefore, Eq. ([T]) implies that, 
for any non-zero value of viscosity, one cannot take the limit of vanishing relaxation time. 
In this sense, the transient dynamics of the dissipative currents cannot be neglected when 
dealing with relativistic fluids. This realization forms the basis for the following discussion. 

In this paper, we derive the equations of motion for dissipative currents in the linear 
regime. We show how the linearized equation of motion of any dissipative current is deter- 
mined once the analytical structure of the associated retarded Green's function is known. 
In the case that the singularity of the retarded Green's function nearest to the origin is a 
simple pole on the imaginary axis, we prove that equation of motion can be reduced to a 
relaxation-type equation. The relaxation time is equal to minus the inverse of the imaginary 
part of this pole. We prove that this prescription gives a value for the relaxation time that, 
under certain simplifying assumptions, coincides with the one derived by matching rela- 
tivistic fluid dynamics to kinetic theory. We also demonstrate that attempts to derive the 
relaxation time from considering the long- wavelength, low-frequency (i.e., fluid-dynamical) 
limit of the retarded Green's function in general fail to give the correct result. This shows 
that transient dynamics is determined by the slowest microscopic and not by the fastest 



fluid-dynamical time scale. 

This paper is organized as follows. In the next section we define the notation used 
throughout this paper. In Sec. IIIII we show the equivalency of the gradient expansion in 
coordinate space with a Taylor series in momentum space. In Sec. IIVI we derive the main 
results of this paper and establish the connection between the analytical structure of retarded 
Green's functions and the linear transport coefficients in fluid dynamics. We show in Sees. 
IV Al and IV Bl that the linear transport coefficients computed using either the linearized 
Boltzmann equation or via the disturbances in the space-time metric follow the general 
method and formulae derived in Sec. IIVI In Sec. |VT] we compare our results to previous 
derivations of relaxation-type equations for the dissipative currents based on the gradient 
expansion. We conclude this paper with a summary of our results. 

II. DEFINITIONS 

Let us consider a general linear relation between a dissipative current J (X) and a ther- 
modynamical force F (X), 



where X = {t, x) is the coordinate four- vector in space-time. Any translationally invariant 
theory that has a linear relation between J and F can always be written in the form of Eq. 
([2]). In fluid dynamics, J could e.g. be the shear stress tensor vr^^ and F the shear tensor 
a'^^, or J could be the diffusion current of a density n and F ~ dn. 
We define our Fourier transformation in the following way 



Here, Q = [u, q) is the momentum four-vector, and Q ■ X = q'^x^ is the scalar four-product 
of the four-momentum vector Q with the coordinate four- vector X. Our metric signature 
is ?7^i/ = diag (+,—,—, —). Using this convention, we can rewrite Eq. (|2]) in terms of the 
Fourier transforms of the retarded Green's function and of the thermodynamic force, Gr (Q) 
and F {Q), respectively, which then gives 




(2) 




(3) 




J{Q) = Gn{Q)F{Q) 



(5) 
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We only consider systems where the microscopic dynamics is invariant under time reversal 
and, thus, ReGj? is an even function of u, while ImGa is an odd function of u. Note that 
Eq. implies that the current J can also be expressed as an integral over Q, 



Thus, Gfi at all frequencies can contribute to the dynamics of J. Equations (E]), (|5]), and 
are, of course, equivalent and contain all the information about the underlying microscopic 
theory that can be obtained through a linear analysis. In this work, we show how the 
analytical structure of the retarded Green's function Gji determines the equation of motion 
for the current J. 

III. EQUIVALENCY BETWEEN GRADIENT EXPANSION AND TAYLOR SE- 
RIES 

In this section, we show that the gradient expansion in space-time is actually equivalent to 
a Taylor series in 4-momentum space. For the sake of simplicity, we suppress any dependence 
on spatial coordinates or, equivalently, on 3-momentum, retaining only the dependence on 
time t and (complex) frequency uj. The coordinate, or 3-momentum dependence, respec- 
tively, will be restored later. We implicitly work in the rest frame of the fluid, such that the 
equations of motion do not appear to be relativistically covariant, however, covariance can 
be restored by a proper Lorentz-boost. 

Let us assume that Gr {u) is analytic in the whole complex u plane. This can be con- 
sidered as the limiting case when Gr (uj) has singularities, but all of them are pushed to 
infinity by some suitable limiting procedure. In this situation a Taylor expansion of Gr {u) 
around the origin has infinite convergence radius and thus provides a valid representation 
of Gr (u) in the whole complex plane. 



Using Eqs. and ([7]), as well as the Fourier representation of Dirac's delta function 






(7) 




(8) 
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it is straightforward to obtain the general form of Gn {t — t') 



Gr {t - t') = Gr (0) 6 {t - t') + id^ Gr {u) dt5 {t - t') 

" ' dl5{t~t') + 0{c^) . (9) 



-^dlGniu) 



'R 

Substituting Eq. OH]) into Eq. we obtain the following equation of motion for the dissi- 
pative current, 

J (t) = 5oF (t) + D^dtF if) + D^dlF (t) + O [dfF] , (10) 
where we introduced the coefficients 

Do = Gr (0) , 



Di = id^ Gr (uj) 



(11) 

uj=0 



D2 = -ldlGRiuj] 



uj=0 

(12) 



This is nothing but the so-called gradient expansion, in which the current J is expressed in 
terms of the thermodynamic force F and its derivatives. Note that the standard gradient 
expansion does not involve time derivatives of the thermodynamic force. This is not a 
problem, since one can always replace time derivatives with spatial gradients using the 
conservation equations of fluid dynamics. 

It is important to remark that when the system exhibits a clear separation between the 
typical microscopic and macroscopic scales, A and i, respectively, it is possible to truncate 
the expansion on the right-hand side of Eq. (ITOl) . A microscopic scale is, for example, the 
mean-free path in dilute gases. A macroscopic scale is given by the inverse of the gradient 
of a macroscopic variable, such as energy density, charge density, or fluid velocity. Note that 
the thermodynamic force, F, is already proportional to a gradient of a macroscopic variable, 
and thus F ~ i^^. Every additional derivative dt brings in another inverse power of £, 
d'^F ~ The microscopic scale A is contained in Gr and its derivatives with respect 

to uo. Thus, up to some overall power of A (which restores the correct scaling dimension), 
Gij(O) ~ A, and each additional derivative dui brings in another power of A, such that 
Dn ~ A"+^ Therefore, the terms DqF, DidtF, and D2dfF in Eq. are of order X/i, 
(A/£)^ and {X/i)^, respectively. This is a series of powers in the so-called Knudsen number 



Kn = A/£. If Kn ^ 1, the gradient expansion of F, Eq. (fTOj) . can be truncated at a given 
order and one obtains a closed macroscopic theory for the dissipative current J. 

IV. THE ROLE OF THE ANALYTICAL STRUCTURE OF Gr {uj) 

In the previous section, it was shown how to relate the gradient expansion of the ther- 
modynamical force F with the Taylor expansion of the Green's function Gr. The viability 
of the latter required the assumption that the singularities of Gr are pushed to infinity by 
some suitable limiting procedure. For instance, if F{uj) has only support in a region of small 
which is well separated from the singularities of Gr, one can devise a limiting procedure 
that effectively pushes these singularities to infinity. However, a priori it is not at all clear 
that F{uj) has vanishing support in the region of the complex uj plane, where Gr {oj) has 
singularities. 

Therefore, we have to consider the case that Gr {uj) has some singularities in the complex 
Lo plane. This fact necessarily restricts the convergence radius of the Taylor expansion. If 
the singularities are simple poles, it is better to use a Laurent expansion around these poles. 

A. Gr {ijj) with one pole 

In order to illustrate this, we consider a retarded Green's function, Gr {uj), with a single 
simple pole at ujq. See Fig. [1] for an illustration in the complex plane. A function with a 
single pole can always be expressed in the following form, 

Gr {uj) = ^ , (13) 

UJ — UJq 

where / {u) is an analytic function in the complex plane. In order for ReG/j to be an even 
function of uj and ImGR to be an odd function of u, we have to require that Uq = 
where ( is positive and real, for the retarded Green's function. We also have to require that 
Re/(w) is odd in u, while lmf{uj) is even in uj. Since J is not a conserved quantity, we 
exclude the case where the pole is at the origin, C = 0. 

Since Gr {uj) has a pole at ujq, the Taylor series around u = has a radius of convergence 
\uo\ and, consequently, this expansion is not able to describe Gr{u) beyond the pole. In 
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such cases, the Laurent expansion for (aj)around the pole uq should be used, 



Gr(uj) 



U! — U!q 



+ f (c^)L=.„ + T;dl f (u;)L=.„ (u; - u^) + O [{u - uof] . (14) 
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FIG. 1. Analytic structure of the retarded Green's function with a singularity in wq- The dashed 
line illustrates the radius of convergence of the Taylor expansion around the origin. 



The series of positive powers in w — can be rearranged into a series of positive powers 
in bj. The coefficients of the latter series can be most conveniently expressed by matching 
the Laurent expansion to the Taylor expansion around w = 0. To this end, we expand 
/ (wq) (w — Wq)^^ around the origin. Then we match Eq. ([7]) with Eq. f lT^ . In this way, we 
obtain e.g. for the coefficients of the constant and linear terms in u 



G«(0) 



Gr [u] 



w=0 



Wo + • • • 



\U)=U)Q 



+ 



(15) 



We now observe that all constant terms in u in Eq. ([H]) can be expressed as Gr{Q) + 
f{coo)/uo, while all linear terms in u can be written as d^jG r{u)\uj=o + fif^o)/^^- Higher- 
order terms in u can be expressed in a similar fashion. Thus, we can rewrite Eq. (IHj) in the 
following way 



Gr{u) 



U — Uq 



+ 



Gr (0) + 



+ 



duj Gr (w) 



aj=0 



+ 



J 



Lo + OM. (16) 



The inverse Fourier transform of Gr (w) is straightforwardly obtained and has the form 



Gr {t - t') = -if M exp [-lUJ^ {t - t')] e {t - t') + 



Gr (0) 



5 it- t') 



+1 



du, Gr (u) 



uj=0 



(17) 



This retarded Green's function is the solution of the following differential equation 



dtGR it - t') + iuj^Gr it - t') = iuoGr (0) 5 (t - t') 



uJoduj Gr (uj) 



uj=0 



Gr{0) dt6{t-t') + 0{d^). (18) 



Dividing Eq. f lTSj) by iuo, multiplying by F {t'), and integrating over t', one now obtains an 
equation of motion for the current J defined in Eq. (|2]), instead of a simple algebraic identity 
as in the gradient expansion, cf. Eq. ( ITOl) . This equation of motion reads 

(19) 



TRdtJ + J = DoF + DAF + O {d^F) 



where the coefficients are 

tr -- 
Do-- 
Di-- 
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iuo C' 
Gr (0), 

id^ Gr (u) 



uj=0 



+ DqTr = TRd^ f (u) 



\uj=0 ■ 



(20) 



Note that, in case Gr{uj) has a single simple pole, Eq. ( [T9|) is exact. 

It is clear that Eq. f|T9|) is nothing but a relaxation equation which is similar in structure to 
that occurring in the transient theories for non-relativistic fluid dynamics proposed by Grad 
10 1 and extended to relativistic fluids by Israel and Stewart in Ref. j^. The appearance 



of the time derivative of J is due to the existence of the pole in the retarded Green's 
function. Also, the transport coefficient tr, usually known as the relaxation time coefficient, 
is directly related to the singularity, ojq = —i(, of Gr{ijj). Since ( > 0, tr is real and 
positive, as expected. It is interesting to note that Dq = Dq, with Dq from Eq. ( ITTj) . while 
Di is not identical to Di, cf. see Eqs. ffTT]) and f l2U]) . 

It should be noted that the right-hand side of Eq. f fTOj) can be truncated using the same 
procedure employed for Eq. (fTOj) . After this truncation, the dissipative current J satisfies 
the following equation, 

TRdtJ + J = DoF + D.dtF . (21) 



As before, it was assumed that Kn ^ 1, with Kn being the appropriate Knudsen number, 
and terms of order O (Kn^) were dropped. 

As was mentioned in the beginning of this section, the Taylor expansion around w = is 
vahd in a radius \ujq\ around the origin, see Fig. [1] Thus, when the pole is pushed to infinity, 
I Wo I — > oo, the radius of convergence of the Taylor series becomes infinite and we should 
recover the gradient expansion. In fact, taking the limit |u;o| — )■ oo, that is, tr — )■ 0, cf. Eqs. 
f l20p . one recovers Eq. fllOl) . with identical coefficients. The coefficient Dq was already seen 



to 



3e identical to Dq, while Di agrees with Di only in the limit of vanishing relaxation time 
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In the case where the thermodynamic force varies slowly on the time scale given by r^j, 
eventually, i.e., for times t ^ r/j, the dissipative current J will follow the time dependence 
imposed by the right-hand side of Eq. ( |2T|) . In other words, the transient term TfidtJ in Eq. 
f l2T]) will become small. Then it is permissible to replace J in this term by the right-hand 
side of Eq. ( 12T|) . and we obtain up to terms of order C(Kn^) 



J-DoF + DidtF - TRDodtF = D^F + D^dtF , (22) 



i.e., we recover the result (1101) given by the gradient expansion. In this sense, the gradient 
expansion is the asymptotic solution of Eq. (fT9l) for time t ^ r/j. 

For nonrelativistic systems, when the transient dynamics can be neglected at all times, 
it is known that the first-order truncation of the gradient expansion can actually serve not 
only as an asymptotic solution, but as an effective theory to describe the system, e.g., sub- 
stituting the first-order result into the conservation equations one obtains the nonrelativistic 
diffusion equation and the nonrelativistic Navier-Stokes equation. However, for relativistic 
theories this is not possible because of the violation of causality. As was mentioned in the 
Introduction, because of Eq. ([1]) for any non-zero value of the transport coefficients (shear 
viscosity, bulk viscosity, etc.) it is not possible to take the (acausal) limit of vanishing re- 
laxation time, if one wants to obtain causal and stable relativistic fluid-dynamical equations 
of motion. 
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B. Gr (oo) with two poles 



In order to better understand the consequences induced by the retarded Green's function's 
nontrivial analytic structure, it is useful to analyze in detail the case where Gr (uj) has two 
poles, uji and 002, 

(23) 



OJ — OJi u — U2 

We employ exactly the same steps as before and expand each term of Gr (uj) in a Laurent 
series around its respective pole. The result is 



U! — U!i U! — U!2 



la;=tJ2 



{CO -iJ2) + [{CO - cOi f] . (24) 



As before, we can match this expansion to the Taylor expansion near the origin. This enables 
us to rewrite Gr {u) as 



G«(u;) = AM + AM + 



CO — COi CO — CO2 



Gr{0) + ^-^ + ^-^ 



C02 



+ 
+ 



Gr {u) 
ld^,GR(oo) 



oj=Q 



h (^1) /2 (^^2) 

,2 , ,2 



COi 



COo 



CO 



uj=0 



fl{(^l) f2{uJ2) 
,3 , ,3 



CO^ 



COn 



co^ + 0{co^). (25) 



With the last expression, we can determine the Green's function GR{t — t') and also its 
equation of motion. Similarly to the previous section, it is straightforward to show that 

Gr {t - t') = -i { /i (a;i) exp [-ico^ {t - t')] + {C02) exp [-ico2 {t - t')] ^ 9 {t - t') 

S{t- t') 



+ 



G^(o) + AM + AM 



Gr iu) 



w=0 



002 



/l(^l) 72(^2) 
,2 , ,2 



cot 



COn 



dlGR{co) 



+ 



/l(c^l) , /2(W2) 



+ 



u)=0 



CO, 



COo 



dtS{t-t') 

dl 6{t - t') + 0{dl) . 



(26) 



The equation satisfied by Gr is 

-dlGR {t -t')^i (a;i + U2) dtGR {t - t') - CO1CO2GR {t - t') + CO1CO2GR (0) S {t - t') 



+ i 



0JlUJ2d^ Gr {u) - {cOr + CO2) Gr (0) dt5 {t - t') 



uj=0 



- ui C02 dl Gr{co) - (a;i + C02) Gr{co) + Gr{Q) 

A a)=0 a;=0 



dl8{t-^)-rO{dl), (27) 



11 



The main difference to the previous case is that, due to the existence of a second pole, the 
Green's function now satisfies a second-order differential equation, instead of a first-order 
one. As will be shown later, the order of the differential equation satisfied by Gji (u) is equal 
to the number of its singularities. Dividing by —U1U2, multiplying the equation by F{t'), 
and integrating over t' we can determine the equation of motion for J, 

X2dU + XidtJ +J = DoF (t) + D^dtF (t) + D^d^.F (t) + O [{X/i)'] . (28) 

We introduced the following transport coefficients, 

1 

X2 = , 

UJ1U2 

1 1 

Xi = -■ — \- 
Do = Gr (0), 

Di = iduj Gr (w) + Doxi , 

D2 = -]-dlGR{uj) +D,xi + D,{x2-xV>. (29) 

Note that X2 and Xi have contributions from both poles. 

Next, we shall investigate under which circumstances a relaxation equation for J can 
be obtained. Due to time reversal invariance, the two poles of Gr (w) can appear in two 
ways: (i) both poles are on the imaginary axis, in which case we assume, without any loss 
of generality, that |ci;2| > (ii) both poles have the same imaginary part, but opposite 
real parts, being symmetric with respect to the imaginary axis. In this case, l^il = \uj2\- See 
Fig. [2] for details. Both cases reflect distinct physical scenarios. 

Let us consider the case in which the thermodynamic force is turned off, and the current 
J is left to relax to equilibrium. This process is governed by the equation 

X2dlJ + XidtJ + J = Q. (30) 

This is the equation of motion of a damped harmonic oscillator, 

x + 2'^x + ulx = Q . (31) 

The harmonic oscillator is overdamped, if 7 > cjq, and underdamped, if 7 < Wq. Identifying 
the coefficients, we find 

7 = ^, ^1 = ^. (32) 

^ X2 X2 
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FIG. 2. Analytic structure of the retarded Green's function with two singularities, uoi and lo2- Case 
(i), in which both poles are on the imaginary axis, is illustrated on the left. Case (ii), in which 
both poles are symmetric around the imaginary axis, is illustrated on the right. The dashed line 
illustrates the radius of convergence of the Taylor expansion around the origin. 

Therefore, if Xi > 4x2? the dissipative current J relaxes to equilibrium without oscillating, 
while for < 4x2, it relaxes in an oscillatory fashion. Using the definitions ( 129|) . we 
see that in the overdamped case, 4:UiU2 > (wi + ^2)^, while in the underdamped case, 
AUJ1UJ2 < (wi + ^2)^- 

In case (i), both poles are purely imaginary, coi = —iQ, with d > 0, i = 1, 2. Since always 
4C1C2 < (Ci + (2)^, we are in the overdamped limit. In this case, a relaxation equation is 
obtained by applying the limiting procedure X2 0. That is, the second pole is pushed to 
infinity, in which case the relaxation time is given by the inverse of the first pole, 

1 

tr = Xi = —■ 

Naturally, this theory is only valid for frequencies that are small compared to the second 
pole. 

In case (ii), l^il = \uj2\ and Ui = — We always have |cjip > (Imwi)^, which can be 
rewritten in the form 4uiU2 = — 4|ci;ip < {ui — w*)^ = {ui + ^2)^. Thus, we are in the 
underdamped limit. In this case, due to the symmetries of the retarded Green's function, 
one cannot disregard one of the poles while keeping the other. Thus, the term including the 
second time derivative of J must be included (otherwise there would be no oscillation) and 
the full equation of motion must be solved. It is important to remark that, in this case, the 
coefficient Xi cannot be interpreted as a relaxation time. 
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Note that the exponential decay of J{t) in case (i) can also be seen when inserting the 
Green's function Gr{u) into Eq. ([6]) and performing the u integral via contour integration, 
picking up the poles via the residue theorem. In case (ii), one gets an exponentially damped 
factor from the imaginary part of the pole, ~ exp(— |Imu;i| t), and an oscillatory factor from 
the real part, ~ exp{iR.eUit). If llmwjl ^ |Rea;j|, the damping is much stronger than the 
oscillation, already within a single oscillation, the dissipative current has decayed a couple of 
e— folds towards its stationary solution. On the other hand, if |Ima;j| -C |Rea;j|, the current 
will oscillate many times before a substantial decay occurs. The oscillatory behavior of the 
dissipative current in case (ii) was already noticed in Ref. 121 ]. 

It is clear from this analysis that the derivation of relaxation equations for systems with 
more than one pole is not possible if the first pole does not lie on the imaginary axis. With 
this in mind, we can finally consider the general case of an arbitrary number of poles. The 
result will be qualitatively similar to what was found above. 



C. Gr (Q) with N poles 

Now we assume that Gr (Q) has poles in the complex w-plane, (q), . . . , Ujy (q). 
We furthermore restore the q dependence which was neglected in the previous sections, and 
assume that the Green's function is analytic in q. Since J is not a conserved quantity, we 
can safely assume that all poles remain within a finite distance from the origin even when 
q — 0. Here we consider the case of a finite (but arbitrarily large) number of poles. Also, 
in many cases Gr (Q) contains branch cuts. We assume that these branch cuts are located 
further from the origin of the complex frequency plane than these poles and thus will not 
be considered any further in this work. The retarded Green's function can be written in the 
general form 

Gr (Q) = — N = 7 r~vi~^ TTi ' (^^) 

j^u-Ui (q) [u - Ui (q)J ■■■[u-Un (q)J 

where the fi (Q) and H (Q) are analytic functions in the complex u plane. We now use the 
identity 

1 + (q) + + (q) = (-1)" l--M^)]-l---.(n)] 

ui (q) ■■■UN (q) 

14 



where 



*m (q) = (-i)' 



E 



1 



(35) 



l<ii...<im<N 

From this expression, we can find the following equation for G^i {Q), 

N 



1 + $1 (q) i-tu) + . . . + <l>jv (q) i-tur Gn (Q) 



-ifEiQ) 



(36) 



Ui [q)--- CON iqj 

After taking the Fourier transform and expanding the functions $m (q) in a Taylor series 
around q = 0, we obtain a differential equation satisfied by Gr. It is clear that this will 
be a linear differential equation of order N in time. The equation of motion for J can be 
obtained in the same way as before. The result is 



XnO^J +... + xidtJ + J = DoF + ... + DmB^F + O {d^+'F, d^) . 

Here, we omitted all the terms involving spatial derivatives. The coefficients Xm are 

1 



(37) 



Xm (0) 



(0) 



(38) 



l<ti...<im<iV 

As before, the coefficients Dq, Di and D2 can be expressed as 







Gr (u;,0) 



Di = idu, Gr {u, 

1 
2 



1^=0 



+ DoX 



1 ; 



D2 = -^dlGR (u,0] 



uj=0 



+ Dixi + Do {x2 - xT) 



(39) 



cf. Eq. fl2^ . In general, the coefficients have the following form 

.(-1)^ dtEico,0) 



\uj=0 



(40) 



k\ wi (0) ■ ■ ■ wtv (0) 

We remark that up to now we have not employed any approximation besides the initial 
assumptions regarding the singularities of Gr. As mentioned before, if there is a clear 
separation of scales, it is possible to simplify the equation of motion. Since Dn ~ A""*"^, 
assuming that X/£ <^ 1, we can safely assume DqF ^ DidtF ^ D2d^F and truncate the 
right-hand side of Eq. fl57|) . 

As discussed in the previous section, a relaxation equation can be obtained only for 
the cases in which the pole nearest to the origin (in the following referred to as the "first 
pole") lies on the imaginary axis. In this case, the limiting procedure which pushes all the 
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other poles to infinity, Xj — > 0, i > 2, can be applied without breaking any symmetries of 
the retarded Green's function. Assuming that this can be done, we obtain the following 
equation of motion for J 



TRdtJ + J^DoF + DidtF + D2d^F + O {Dsd^F, d^) . 



(41) 



where 



1 



iui (0) ' 
Do = Gr {co, 0) 

uj=0 

Di = id^ Gr (u, 0) 

D2^-ld'^GR{u,0) 



uj=0 



LJ = 



DqTr , 

+ DiTr-DotI , 



(42) 



On the other hand, if the first pole and, for reasons of symmetry, its counterpart on the 
other side of the imaginary axis, have nonzero real parts, the dissipative current will oscillate 
and the equation of motion cannot be reduced to a simple relaxation equation, even in the 
small-frequency domain. 



V. APPLICATIONS 



A. The Linearized Boltzmann Equation 

The discussion presented above is valid for both weakly and strongly coupled theories. 
The prime example of a weakly coupled theory is given by the Boltzmann equation. In this 
section we calculate the shear viscosity and relaxation time coefficients for a weakly coupled 
gas via the Boltzmann equation following the general method presented above. 

We start from the relativistic Boltzmann equation 

K-dU^G [/] , (43) 

where K = (/co,k), ko = + m?, and m is the particle mass. We use the notation 
/k {^) — f {X, K) . We consider the linearized Boltzmann equation around the classical 
equilibrium state, /ok = exp (yok), where yok — ocq — PoEyg^, with the inverse temperature Pq, 
the ratio of the chemical potential to temperature ckq, and the energy in the local rest frame 
E\^ = U ■ K, respectively. Here, is the fluid four-velocity. 
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The Boltzmann equation can be written as 

6h + E^'K ■ V6h - C {X, K) 6h = S (X, K) , (44) 

where we defined 5/k = /k — /ok, A = U ■ dA, = A^^d", with A^" = r]^" — u^u"^ being 
the 3-space projector orthogonal to m^, and 

S (X, K) = /ok [-ao + Ejo + K ■ [fi^U + V/3o - E^^Va, 



(45) 



2 — K V""" — K/ - ' i-^ — K 

Here, we introduced the expansion scalar 9 = d ■ U and the shear tensor a^'^ = V^^u'^\ 
where = A^"'''^A^f} with A^"'''^ = (A^^A'^^ + A'^^A^") /2 - A^'^A^^/S being the 

double symmetric traceless projection operator. In Eq. (14^ . we also introduced the collision 
operator 

C (X, K) Sh = ^ f dK'dPdP'W^^,^^^Jo^.foi.' f ^ + ^ - ^ - ^"j , (46) 

'^-c^k J V /Op' JOp JOk' JOk / 

where VTkk/^pp/ is the transition rate, u is the symmetry factor (z/ = 2 for identical particles) 
and 

rfi^ = . (47) 

We now consider Eq. (jH]) in the local rest frame, where U = (1,0,0,0). We assume a 
situation where the fiuid does not accelerate, [7 = 0, and does not expand, = 0, and where 
temperature and chemical potential are constant. With these assumptions the source term 
reduces to 

S (X, K) = f3oUE^'k^''k''^cr,, (X) , (48) 

and the collision operator no longer depends on X, C{X,K) = C[K). The Boltzmann 
equation (jUj) takes the form 

9t5/k + V ■ V<5/k - C (K) 6h = S (X, K) , (49) 

where v = h/E]^. 

We solve the inhomogeneous, linear, partial integro-differential equation (H9ll for (5/k in 
four-momentum space. Taking the Fourier transform we obtain 

- too6f\ iQ) + w q 5/k (Q) - C {K) 6f\ {Q)=S{Q,K) . (50) 
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From now on, it is important not to confuse the momenta of the particles, K, K', P, and 
P', with the variable Q from the Fourier transformation. 

The formal solution of Eq. (1501) can be expressed in the following form 

iQ) = — — (Q, K) . (51) 

—lu + zv • q — C (A ) 
Note that the shear stress tensor, tt^", can be expressed in terms of 5/k as [l^ 

dKk^^'k''^5h. (52) 

Taking the Fourier transform of Eq. (!52l) and using Eqs. (HHj) and (ISTj) . we obtain 

(Q) = / dRk^'^k"^ i B^f^^E^^k^-k^)a^^ (Q) . (53) 

J —iu + zv ■ q — C7 [K) 

Note that Eq. (153|) has the following form, 

fr'^- (Q) = Gr'' (Q) (Q) , (54) 

where we introduced 

GT^ iQ) = [ dKk^^k'^) ^ ^^^i5^f^^E-'k^^k^\ (55) 

J —loj + ■ q — C (a ) 

As already mentioned, we want to calculate the two main transport coefficients that 
describe the linearized dynamics of the shear stress tensor: the shear relaxation time, r^r, 
and the shear viscosity coefficient, rj. Let us start by defining the function 

B^P = — ^^^p,f,^E^^k^-k^\ (56) 

—luj + «v ■ q — C (A ) 

where, by definition, 3°"^ satisfies 

-lu + iwci-C (ir)] B'^^ (g, K) = PoE^'k^'^k^^fo^. (57) 

In general, B""^ is a function of Q and K. However, from Eqs. fj38|) and fl39|) we already know 
that, in order to calculate the relaxation time and the viscosity coefficient, it is sufficient to 
consider the case q = 0. Then = {uj,K). The dependence of on K can be 
expressed via the following expansion, 

oo 

B^^ {oj, K) = Uk^^k^^ Yl (^) ^k- (58) 



n=0 
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Substituting Eq. ( 158|) into Eq. (155|) . it follows that 

oo „ 



n=0 

oo 



2A'^'^-^^J„+4,2anM , (59) 



n=0 

where we used LK 



y dKk^^k''h^''k^'>E'^U^ = ^^ A^^-^^ j dKE^U^{m^-Elf , (60) 
and introduced the thermodynamic function 

In, = I dKf,^E-'^ (m^ ~Eiy. (61) 

Thus, Eq. ( I5^ can be cast into a more convenient form 

n^'" {oo, 0) = 2Gr {oo, 0) a^"" {oo, 0) , (62) 
where we introduced the retarded Green's function 

oo 

(63) 

n=0 

Since in fluid dynamics the thermodynamic force related to the shear stress tensor is defined 
as 2a^'^ , we kept the factor 2 in Eq. (|62l) . As shown in the previous section, Eqs. (l38l) and 
( l39l) . the shear relaxation time is determined by the first pole, ooi, of Gr {oo, 0), as 

= (64) 

looi (0) 

while the shear viscosity coefficient is determined by the retarded Green's function Gr {oo, 0) 
at the origin, 

7]=Gr {oo, 0) . (65) 

Thus, the problem of finding the linear transport coefficients has been reduced to determining 
the analytic properties of a„ {oo). 

Equations such as Eq. (!57|) appear quite often in problems involving the extraction of 
transport coefficients from the Boltzmann equation. The way to solve this problem is to 
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substitute the expansion fl^ into Eq. f lFTI) . multiply by E^k^^^k'^\ and integrate over dK. 
Then one obtains 

oo „ 

^a„(u;) I dKE^k'^^k"'^ -loo - C (K) foi^E^^k^'' k^^ 

= (3oj dKEl^-^k^^^k'^k^^k^^foi, , (66) 
Using Eqs. (1601) and (16T|) . we can rewrite this in the form 

oo 

J2 (-'^cul?'™ + A'^n a„ (a;) = /3o/™+3,2 , (67) 



n=0 

where we defined the matrices 



A'^'^A^""'^ = -\j dKE^k'^^k^'^C (K) fokEl^k'^'' k^\ (68) 

= ^ / ^^^01^^^^" ("'^ - ^k)' • (69) 
Thus, the formal solution for an{uj) is 

oo 

am (w) = /3o 5Z [i-tuV + A)-'] 4+3,2 (70) 

n=0 

and the expression for Gr {u, 0) becomes 

oo oo 

Gn {oo, 0) = PoY,Yl ^-+4.2 [i-'^'^ + ^)"'] (71) 
m=0 n=0 

The poles of the function above can be obtained from the roots of the determinant 

det {-icoV + ^) = , (72) 

which is given by the product of the eigenvalues A„ of the operator {—iuV + A). Note 
that, because V and A are real matrices, all the poles are on the imaginary axis. Thus, the 
truncation of the equation of motion to a relaxation-type form is possible, if the separation 
between the poles is large enough. The shear viscosity coefficient is given by 



oo oo 



V^Gr{u,0) _ =PoJ2Yl W4,2(^-')"^"/n+3,2 ■ (73) 



a;=0 

m=0 n=0 



Thus, in order to find the relaxation times and viscosity coefficients from the linearized 
Boltzmann equation one has to invert and compute eigenvalues of infinite matrices. In 
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practice, however, one never deals with infinite matrices because the expansion ( !58|) is always 



truncated, see, e.g., Chapman-Enskog theory 



13|- 



Let us consider the simplest possible case and take only one term in the expansion 
We remark that this corresponds to using the Israel- Stewart 14-moment approximation in 
the moments method. Then, Gr (cj, 0) has the following simple form 

where we used that = I^i- In this case, the retarded Green's function has only one pole, 
uq = —iJ^^/Ii2- The relaxation time is obtained as 



On the other hand, the shear viscosity is given by 77 = (w, 0) 



and becomes 



In the massless limit, for a gas of hard spheres, one determines A^^ to be 

A''^ = ^h2noa, (77) 
5 

where a is the total cross section and no is the particle number density. Then, 

' (79) 



where we used that, in the massless limit, 

^42 = 4§ , (80) 
Po 

^32 = , (81) 

where Pq is the thermodynamic pressure. Also, in this particular example, one can show 
that the ratio r^/r^r is independent of the cross section, 

- = /9o/32 . (82) 

These are exactly the results obtained in Ref. [l^. This demonstrates that the relaxation 
time in IS theories, which determines the time scale of the transient dynamics of the dissipa- 
tive currents, is indeed a microscopic and not a fluid-dynamical time scale. It is determined 
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by the interparticle scattering rate, and not by the time scales of fluid dynamics located 
near the origin of the complex cj— plane. We shall demonstrate in Sec. |Vl]that attempts to 
extract the value of the relaxation time from the dynamics on fluid-dynamical time scales 
in general fail to give the correct expression. 

It is important to remark that by including more terms in the expansion (135]) we obtain 
a retarded Green's function with more poles. Furthermore, the expression for the first 
pole and, consequently, the relaxation time, will also be modified. All the other transport 
coefficients will also receive corrections. 



B. Linear Response Theory and Metric Perturbations 



We now apply our formalism to the case studied in Ref. 15|, where the transport coeffi- 



cients are determined from perturbations h^^" of the metric tensor 

gf^^ = n^^^^ + h^^-_ (83) 

This method can be equally applied at strong and weak coupling.^ The variation of the 
energy-momentum tensor T'^'^ due to the metric perturbations is [16] 

ST^" (X) = - / d^X' G^j^""^ {X - X') {X') , (84) 

2 J ~oo 



where G^""^ (X — X') is the retarded Green's function. 

For the sake of simplicity, a very peculiar type of metric perturbation is considered, 
hxy = hxy{t,z), with all other components of the metric tensor left unperturbed jl5, [l7 |. 
For this specific type of metric perturbation, all the other components of ST'^'^ decouple from 
the xy component and one arrives at a very simple expression for ST^^ 1^ 

/oo 
dt' dz' G^^"^ (t -t']z- z') hxy {t', z') . (85) 
-oo 

The energy-momentum tensor T^'^ is then assumed to have the traditional fluid-dynamical 
structure 

T^"^ = e - A^'T + -Kf"" . (86) 

For the particular metric perturbation considered here, the effects of bulk viscosity can be 
neglected without loss of generality. It is possible to show that 
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where Stt^^^ is the xy component of the shear stress tensor generated by the metric pertur- 
bations and Pq is the pressure of the unperturbed state. For the sake of simphcity, the shear 
stress tensor of the unperturbed state is set to zero. Also, it is easy to see, using the equa- 
tions of motion, that the velocity terms in the energy-momentum tensor do not contribute 
to 6T^y. Thus, we arrive at the following equation 

/oo 
dt'dz'GY'' (t-t'; z-z') h,y {t' 
-oo 

After taking the Fourier transform we can show that 



(88) 



5^^y{Q) = Gr{Q) Ky{Q) , (89) 

where Gr{Q) = —Pq + G'^j^^^iQ). Note that h^^ = —h^y. Since the pressure has no 
dependence on Q, it is clear that Gr{Q) has the same analytic structure as G^j^^^iQ). 

Assuming that G^^^^ has poles, Ci;j(q), and that the first pole is purely imaginary, we 
can apply Eq. (1411) and obtain the equation of motion for (Jvr^^, 



T^dtdTl'^y + (Jvr^^ = DoKy + DMxy + D2dlKy + O {dfh^y, d%y) . 



(90) 



Note that if the first pole is not purely imaginary, a simple relaxation equation would not 

be able to describe the transient dynamics, even for long times. In Eq. fl90|) . we introduced 

the following transport coefficients, 
1 



" zui (0) ' 
Do = Gr (cu, 0) 



w=0 



-Po+GT'{^,0) 



= -Po + Po = 



Di = td^ Gr (a;, 0) +t^Do = td^ Gr (cu, 0) 

uj=0 



uj=0 



V 



Do 



1 

—i 
2 



-^dlGR(Lo,0] 



+ DiT^ - DotI 



(91) 



1 

The expression for the shear viscosity rj is now different from the one in the Boltzmann case, 
because 2cr^^ = dth^y, so that the coefficient of the shear tensor is actually Di (and not Dq, 
as before). On the other hand, the relaxation time is still given by the inverse of the first 
pole. 



VI. DISCUSSION 



Relaxation-type equations for the dissipative currents have been recently derived in Ref. 
1^ . We give a brief account of the strategy employed in that work in terms of our notation. 
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The starting point is the gradient expansion ( ITO!) . assuming that a Knudsen number counting 
as explained at the end of Sec. [TTll is apphcable. The gradient expansion ffTOj) is not an 
equation of motion for the dissipative current, but one can construct one by taking the 
first-order solution, J = DqF + 0{Kn^), and then replacing the first time derivative of F 
on the right-hand side by a time derivative of J, 

J = DoF + Didt(^-^^ +0{Kn') ^ fRdtJ+J-DoF, (92) 

where ff> = —Di/Dq has the dimension of time. By construction, this is a relaxation-type 
equation of motion for the dissipative current J, with a relaxation time tr. 

From our previous discussion, however, it is clear that this need not be the correct 
equation of motion for the dissipative current. If the poles of the retarded Green's function 
Gr{u) associated with the dissipative current J are off the imaginary axis, the equation of 
motion for J is never of relaxation type, and the above way to construct one is misleading, 
since it fails to capture the correct physics. Only if the first pole of lies on the 

imaginary axis and is sufficiently separated from the other singularities, one can obtain a 
relaxation-type equation for J. In this case, however, the true relaxation time is tr = 
l/[iui{0)], and not tr = —Di/Dq. 

There is, however, a particular case, where tr = tr, namely when Di = 0, i.e., when 
S(u;,0) = const.. This is most easily seen in the one-pole case, where, cf. Eq. (1201) . 



= td^GR{oo,0) 



+ DoTr = D^ + Do tr , (93) 

u}=0 



i.e., Tr = -Di/Do = Tr. 

This argument can also be applied to the case discussed in Sec. IV B| where the thermody- 
namic force is not given by F but by dtF. Then, Dq = 0, and Di = Di = rj is the transport 
coefficient. In this —D2/D1, and equivalency to the true relaxation time requires 

that D2 = 0, i.e., 

+ D,T^ = D2 + D^T^ , (94) 



= -ldlGR{u,O] 



uj=0 



or 

1.2 



2 



Vr^ = 7;dtGRiuj,0) . (95) 



u}=0 



This equation is rather similar to the one given in Refs. [15|, ll7|, tjTt, = [d'^GR^O] — 



dn^GR{0)]/2 18|. Here, the additional derivatives with respect to momentum enter be- 



cause second-order time derivatives also arise from space-time curvature in Eq. f fTOj) (they 
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were not explicitly denoted in that equation), cf. Ref. [l5|. These are not subjected to the 
construction of a relaxation-type equation for J as explained above. In turn, r/fv receives 
an additional contribution from second-order spatial derivatives, for details, see Ref. 15|. 

Regardless of these considerations, the true relaxation time is always given by the first 
pole of the retarded Green's function. In general, the location of this pole cannot be found 
from a truncated Taylor expansion around the origin. 

A calculation of the shear viscosity relaxation time coefficient has also been performed in 

" n 



Refs. 



19| within the Mori-Zwanzig formalism ^] . The implications of our findings to those 



works will be discussed in detail elsewhere. 



VII. SUMMARY 

In this work, we have derived equations of motion for the dissipative currents, assuming 
these currents to be linearly related to the thermodynamic forces. We have shown how these 
equations of motion are determined by the analytical structure of the associated retarded 
Green's function in the complex u plane. We have demonstrated that the standard gradient 
expansion is equivalent to a Taylor expansion of the retarded Green's function around the 
origin in the complex u plane. This Taylor series is convergent only when all singularities of 
the retarded Green's function are pushed to infinity. In general, however, these singularities 
appear at finite values of |ci;|, which consequently severely restricts the applicability of the 
gradient expansion. 

We have furthermore demonstrated that, if the retarded Green's function has simple poles 
in the complex u plane, the dissipative current obeys a differential equation with source 
terms which are the thermodynamic force and gradients thereof. This is different from 
the gradient expansion where the current is directly proportional to the thermodynamic 
force and its gradients. In general, the equation of motion for the dissipative current is 
not a relaxation-type equation. However, in the limit where all singularities of the retarded 
Green's function except the pole nearest to the origin are pushed to infinity, it is possible 
to approximate the dynamical equation satisfied by J as a relaxation-type equation, similar 
to the ones appearing in Israel- Stewart theory. This is only possible if the first pole is 
purely imaginary. The relaxation time is equal to minus the inverse of the imaginary part 
of the pole. The gradient expansion constitutes the asymptotic solution of the resulting 
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relaxation-type equation, and can be obtained by taking the relaxation time to zero or, 
equivalently, pushing the first pole to infinity. This is consistent with the above statement 
that the gradient expansion arises from a Taylor expansion of the retarded Green's function 
around the origin. 

In relativistic systems, in order to have stable and causal equations of motion for the 
dissipative currents one cannot take the relaxation time to be arbitrarily small [sj. Thus, 
one cannot push the first pole of the retarded Green's function to infinity. This is the reason 
why one cannot use the gradient expansion to obtain stable and causal equations of motion 
for the dissipative currents. 

As an example, we have studied the Boltzmann equation for a classical gas and demon- 
strated that the retarded Green's function associated with the shear stress tensor has in- 
finitely many simple poles along the imaginary axis. Therefore, it is possible to reduce the 
equation of motion for the shear stress tensor to a relaxation-type equation. Our results 



are consistent with those of Ref. 



14l |. when one truncates the collision operator at lowest 



order. This convincingly demonstrates that the time scale of transient dynamics determined 
by the relaxation time is of microscopic, and not of fiuid-dynamical origin: it is the slowest 
microscopic time scale, not the fastest fiuid-dynamical time scale as implicitly assumed in 
attempts to extract the relaxation time by expanding the retarded Green's function around 



the origin. Consequently, the expression for the relaxation time derived in Refs. [15|, [ITj] by 
using an expansion around the origin is, in general, different from the one derived from the 
first pole of the retarded Green's function. In fact, as we have demonstrated in this work, 
when the retarded Green's function has simple poles off the imaginary axis in the complex 
u plane, the true dynamics of the system at long wavelengths and low frequencies is not 
even of relaxation type. Strongly coupled theories, like those emerging from the AdS/CFT 
correspondence |l6|, belong to this class. 
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